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Abstract 



As far as entanglement is concerned, two density matrices of n particles 
are equivalent if they are on the same orbit of the group of local unitary 
transformations, U(di) x ■ ■ ■ x U (d n ) (where the Hilbert space of particle r has 
dimension d r ). We show that for n greater than or equal to two, the number 
of independent parameters needed to specify an n-particle density matrix up 
to equivalence is H r d% — Y2 r d^, + n — l. For n spin-i particles we also show how 
to characterise generic orbits, both by giving an explicit parametrisation of 
the orbits and by finding a finite set of polynomial invariants which separate 
the orbits. 
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In this paper we take some further steps towards understanding multi-particle entangle- 
ment by analysing the non-local properties of density matrices of n particles. This continues 
the programme, begun in |1|], in which we gave a framework for studying the space of pure 
states of n spin-1/2 particles. As discussed in the space of pure states of n spin- 1/2 
particles is (D 2 ™ = C 2 <8> ••• ®C 2 ; however not all the 2 n complex parameters have non-local 
significance: the group of local transformations, U(2) n acts on the space of states and two 
states which may be reached from each other by local actions are equivalent as far as their 
non-local properties are concerned. Each equivalence class of locally equivalent density ma- 
trices is an orbit of this group. For many purposes, only parameters describing non-local 
properties are significant; an example is that any good measure of entanglement must be 
invariant under local transformations, and thus it should be a function of non-local param- 
eters only (here and henceforth we will refer to parameters which are invariant under local 
transformations as invariants). A key question is to identify the invariants. 

In this paper we will focus on density matrices and show that for n > 2, of the 2 2n — 1 
real parameters describing density matrices of n spin-1/2 particles 2 2n — 3n — 1 are invariant 
under local transformations, U(2) n . This generalises to an arbitrary set of n particles as 
Fir $ — J2 r d% + n — 1 where d r is the dimension of the state space of the rth particle. For n 
spin- 1 particles we also show how to characterise generic orbits, both by giving an explicit 
parametrisation of the orbits and by finding a finite set of polynomial invariants which 
separate the orbits. Thus given two density matrices we can compute explicitly whether 
they are on the same orbit or not. Other authors have also discussed the use of invariants 
in discussing entanglement @ || and applied invariant theory to quantum codes ||]. 

In order to calculate the number of functionally independent invariants it will be con- 
venient to find the dimension of the orbit of a generic density matrix under the group of 
local transformations. The dimension of the orbit is the number of parameters describing 
the location of a density matrix on the orbit. The total number of parameters (2 2n — 1 
real parameters) describing the space of density matrices minus the number of parame- 
ters describing a generic orbit (the dimension of the orbit) gives the number of parameters 
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describing the location of the orbit in the space of orbits, i.e. the number of parameters 
describing the non-local properties of the density matrices. 

To fix notation, it will be convenient to consider the case of a one-particle density matrix 
first. The space of pure states of a single spin-1/2 particle is(D 2 and thus a density matrix is 
a 2 x 2 complex matrix which is hermitian, positive and with trace one, and may therefore 
be described by three real parameters. A particularly convenient representation of such a 
matrix is 



P = -J-2 + OLiGi 



where a,-, i = 1, 2, 3 are real and 
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We note that 



V 01 / 



(i) 



(2) 



OiOj — ieijk&k + Sij. 



(3) 



Under a local transformation by a unitary matrix U, p is transformed as 



p i-> UpU ] . 



(4) 



The group £/(2) is isomorphic to U(l) x SU(2) where, physically, the U(l) is the phase 
transformation of a state, represented by a unitary matrix e 1 *^. This element clearly leaves 
any density matrix invariant under the transformation (Q) so that when considering the 
action (^) we may restrict attention to elements of SU(2). In order to find the number of 
invariants it will be more convenient to find the dimension of a generic orbit under the action 
of SU(2). To do so one may work infinitesimally. Thus, associated to the action of the Lie 
algebra of the group of local transformations acting on the space of density matrices there is 
a vector field: if we take an element T of a basis for the Lie algebra, the action of the group 
element k = exp ieT e K on an element p induces an action on functions from p to (D; and 
the vector field, Xt, associated to the Lie algebra element T is found by differentiating: 



XTf(v) d ^§- e f(e- T p)U = §- e f(p + Sp)\ e=0 . (5) 

The linear span of vector fields at the point p associated with the whole Lie algebra forms 
the tangent space to the orbit at the point p and so the number of linearly independent 
vector fields at this point gives the dimension of the orbit. 

A general element of the Lie algebra in the spin- 1/2 representation is given by 

T = r}iOi (6) 

and its action on the density matrix is to give an infinitesimal transformation 

8p = iF,p] (7) 

where [, ] is the matrix commutator. 

We may therefore calculate the three vector fields Xi,X 2 and X 3 associated to the Lie 
algebra elements a%, a 2 and cr 3 as 

d d d d d d 

Xi = «2 -k "3-^ — , X 2 = a 3 - ai- — , A 3 = ofi- a 2 - — • (8) 

oaz oa 2 oa\ oa^ oa 2 oa\ 

We note that at generic values of ai,a 2 ,a3 only two of these vector fields are linearly 
independent since 

aiXi + a 2 X 2 + a 3 X 3 = 0. (9) 

Thus the dimension of the generic orbit is two and therefore of the three parameters describ- 
ing a generic density matrix, two are non-invariant leaving only one invariant parameter, 
as one expects since only the single independent eigenvalue of p is invariant under local 
transformations. 

We note that the effect of the transformations (|4]) is to act on the vector a by rotation 
by an orthogonal matrix, i.e. an element of 50(3) - this follows from the fact that c^ai 
is the representative of a Lie algebra element and the conjugation action (^) is the adjoint 
action of the group on its Lie algebra. We may thus find a way of exhibiting the invariant 
under local transformations: 



I = aiCijdij = \a\ 2 (10) 

where we have used the fact that SO (3) has an invariant tensor 5 tj . We note that this 
invariant may also be expressed as 

I = Tr(p 2 ) - \ (11) 

We now turn to the case of two-particle density matrices. Such a density matrix has 15 real 
parameters, and the maximum dimension that a generic orbit could have is 6 (corresponding 
to two copies of 577(2)) if all the vector fields corresponding to a basis of the Lie algebra 
were independent. We will show that the vector fields do indeed span 6 dimensions, and 
thus that there are 9 non-local parameters. 
We may write a density matrix as 

p = ^12 <S> h + OLi^i <S> h + Pih <S> Oi + RijUi ® Oy (12) 

The action of a Lie algebra element of the subgroup SU (2) acting on the first component of 
the tensor product is 

(jWp = i[r] k cr k <g> l 2 ,p] 

= OL k ri m e mki Oi ® 1 2 + RkjVmtmkiO-i ® CTj, (13) 

and that corresponding to a Lie algebra element of the subgroup SU (2) acting on the second 
component of the tensor product, 

5 (2) p = i[r]kh ®cr k ,p] 

= PkVmtmkih ®CTi + RikVmtmkjVi ® <Tj. (14) 

The vector fields corresponding to the six basis elements <7j <g> I2, I2 ® &i are 



Consider the set Xi first: one can see that these three are linearly independent at generic 
points by considering the coefficients of d/da iy since a linear relation would have to be of the 
form akXk = 0, but one can see that this relation will not hold for non-zero a's by looking 
at the coefficients of the partial derivatives with respect to Rij. Similarly by considering the 
coefficients of the partial derivatives with respect to (3i, (3 2 , (3 3 , one sees that Y 1: Y 2 ,Y 3 are 
linearly independent. Finally, we note that the coefficients of the partial derivatives with 
respect to (3±, (3 2 , (3 3 are zero for X 1 ,X 2 ,X 3 and the coefficients of the partial derivatives 
with respect to a±, a 2 ,a 3 are zero for Yi, Y 2 ,Y 3 so that there can be no linear relation at 
all between the six vector fields X±, X 2 , X 3 , Y\, Y 2 , Y 3 . Thus the dimension of the orbit of a 
generic density matrix is 6 and thus the number of non-local parameters, 15 — 6 = 9. 

In general, we can consider a system of n particles with individual state spaces of di- 
mensions di, . . . , d n . The density matrix is a hermitian D x D matrix with trace 1, where 
D = d\d 2 . . .d n , and therefore requires D 2 — 1 real parameters which can be taken to be the 
coefficients . . . , a^ n \ . . . , R in an expansion 

P = ^U®--"l^ + E«l r) l®---^ (16) 

r=l 

where (i r — 1, . . . , d 2 . — 1) are a basis set of traceless hermitian d r xd r matrices (generators 
of SU(d r )). The action of an infinitesimal generator of SU(d r ) acting on the rth factor of 
the tensor product is 

S (r) P = 4W(4 h ® • • • ® T i ® • • • ® 1 + • • •) k = 1, ■ ■ ■ , d 2 r - 1) (17) 

where c[^ k are the structure constants of SU(d r ). Thus the infinitesimal action of local 
transformations is given by a set of vector fields 

Similar considerations to those used above for the case of two spin-1 particles show that 
these vector fields are generically all independent. Thus the generic orbit has dimension 
^1 "I — ■ d\ — n. Since the space of density matrices has dimension d\ . . . d^, there are a total 
of 
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Y[d 2 r -J2d 2 r + n-l (19) 

r r 

non-local invariants. 

Let us now return to the case of n > 2 spin-1/2 particles and explicitly identify a set of 
invariant parameters which characterise generic orbits. To be explicit, consider the case of 
three spin-1/2 particles with density matrix which may be written as 

P = <S> h <S> h + oiiOi <S> h <S> h + Pih <S> <Ji ® I2 + lih ® I2 ® 

o 

+ RijOi <g> Gj <g> 1 2 + SijOi <g) 1 2 <g> Oj + Tijl 2 ® (7j ® CTj 

+ <5ijfcO"i ® <7j <g) <7fc. (20) 

The action by a local unitary transformation on the first component in the tensor product 
induces the following transformations on the components of p 

oil 1 > LijQLj] Rij 1 ► LikRkj] Sij 1 > LikSkj] Qijk 1 > Li m Q m j k (21) 

where Ljj is an orthogonal matrix, and the other components of p do not change. Similarly 
actions by a local transformations on the second and third components of the tensor product 
induce 

Pi 1— > MijPj] 1— > MjkRik, Tij 1— > MikT k j; Qij k 1— > Mj m Q irnk (22) 

and 

7i l— ^ijlj'i $ij l— *■ Nj k S ik ; 1— > Nj k T ik ; Q ijk 1— > N km Qij rn (23) 

respectively, where and iV^- are orthogonal matrices independent of L. 

We max fix a canonical point on a generic orbit as follows: firstly let us define 

-^Qi' QijkQi'jkj Yjj' QijkQij'ki Z kk i QijkQijk' j (^^) 

and perform unitary transformations on particles 1, 2 and 3 so as to move to a point on the 
orbit in which X, F and Z are diagonal; generically the diagonal entries are distinct and we 



can arrange them in decreasing order (X, Y and Z are hermitian, positive matrices). The 
only remaining transformations which leave X, Y and Z in these forms are local unitary 
transformations which induce orthogonal transformations in which Lij, Mij and JVjj are one 
of the matrices 

\ 



^10 '* 



0-10 
0-1 



^ -1 \ ( -I 



\ 



J 



\ 




0-10 
1 



(25) 



1 
0-1 

We may specify a canonical point on the generic orbit uniquely by specifying that all the 
components of a have the same sign, and similarly for j3 and 7. This method works as long 
as X, Y and Z have distinct eigenvalues and the components of a, f3 and 7 are not zero 
at the canonical point on the orbit. The parameters which describe the generic orbits are 
the components of a, (3, 7, R, S, T and Q at the canonical point on the orbit. We note that 
the number of parameters describing the canonical point are the 2 6 — 1 = 63 components 
of a, P, 7, R, S, T and Q minus the 3x3 = 9 constraints that the non-diagonal elements of 
X, Y and Z are zero; thus the number of non-local parameters is 54 as given by the general 
formula. 

We note that the fact that the canonical point, as constructed, is unique means that 
all points on the same orbit will have the same canonical representative: conversely, if two 
density matrices p\ and p 2 have the same canonical form, then 



UlPlUl = Pcanonical = U 2 p2^ 



(26) 



for some U 1 and U 2 , so that 



P2 = {uiu^pMuiV 



(27) 



and thus p\ and p 2 are on the same orbit. 

We now describe a finite set of polynomial invariants which separate generic orbits by 
finding a set which allows one to calculate the components of a, f3, 7, R, S, T and Q at this 
canonical point. The complete infinite set of polynomial invariants is found by contracting 



the indices of a, /3, 7, R, S,T and Q with the invariant tensors 5ij and e^. However we 
may find a finite set of invariants which separates generic orbits. Firstly we note that 
tr(JT), tr(X 2 ) and tr(X 3 ) determine the diagonal elements X\ , X 2 and A§ of X, and similarly 
for Y and Z. Now consider the three invariants A 2n = a T X n ~ 1 a, n = 1,2, 3. We may write 
these three invariants in the following way: 



\ 



\4j 



U 2 \ 



A, 



(28) 



1 1 1 

x\ xi xl 

\ A2 A3 J 

where ai, a 2 and a 3 are the components of a at the canonical point on the orbit. The 
Vandermonde matrix 



A 



I X 1 ^ 



X\ X\ A3 

i4 \4 \4 



(29) 



V A x A 2 A 3 J 

has determinant {X\ — A2XA2 — A3) (A3 — Xf), and we may solve for a\, a 2 and a\ as long as 
det A is non-zero. Also if the invariant 



Ag = eijkai(Xa)j(X 2 a)k = a.(Xa) A (X 2 a) = aia 2 a 3 det A 



(30) 



is non-zero, then we may determine the sign of the components of a; recall that, by definition, 
all the components of a have the same sign at the canonical point. The analogous expressions 
B 9 , Cg determine the values of (3 and 7 at the canonical point. The values of the components 
of R at the canonical point may be calculated from the following nine invariants: 



I r , s = (X^aUY^P)^. 



These nine equations may be put together into a matrix form 



(31) 



/ = ((AF) ® (MG)) R, 



(32) 



where / and R are column vectors with nine components and the matrices A, F, M and G 
are 
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A 



t y i i\ 



\ A 2 A 3 

\^ ^1 ^2 ^3 J 



( Oi ^ 

a 2 

y o o a 3 y 



M 



A*i A*l 

4 4 4 

\ a*i /4 y 



G 



^ 6i ^ 

b 2 

y o o 6 3 y 



(33) 



where n\, //| and //| are the diagonal elements of Y. We note that det(AF) = A 9 and 
det (MG) = B 9 , so since we are assuming that these are non-zero we may invert the matrix 
equation to find the components i?^. The components of S and T may be found in a similar 
way. Finally we may use the 27 invariants 

Ir,s,t = (X^aUY^P)^-^)^. (34) 

to find the components of Q at the canonical point on the orbit in terms of the I ryS j (there 
will, of course, be some relations between these components due the constraints that X, Y 
and Z are diagonal). 

Thus, by showing that the following set of polynomial invariants is sufficient to calculate 
the components of a generic density matrix at the canonical point we have demonstrated 
that they characterise generic orbits: 

trX r , trF r , trZ r 

a T X r - 1 a, (FY'^P, ^Z^r, 

a.(Xa)A(X*a), @.(Y(3) A (Y 2 [3), 7 .(Z 7 ) A (Z 2 7 ) 

(X^aUY*- 1 ?)^, (Y^PUZ*- 1 ^, (X^aUZ^)^ 

(X»- 1 a) i (y- 1 /3) J -(Z t - 1 7)fcQ<ifc; (35) 

the indices r, s, t range over the values 1, 2, 3. 

If two density matrices have different values of any of these invariants they are not on 
the same orbit; if they have same value of all of these invariants, and if Ag, B 9 and Cg are 
non-zero, then the density matrices are locally equivalent. 

We note that the number of independent components of a generic density matrix at the 
canonical point is equal to the nubmer of functionally independent parameters calculated 
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at the beginning of this letter. However, the number of polynomial invariants needed to 
characterise the generic orbit is greater than this; this is related to the fact that the ring of 
invariants is non-polynomial, i.e. that the geometry of the space of orbits is non-trivial. 

The procedure given above can be used for all n > 2: use the tensors of highest rank 
and rank one in the expression for p to fix a canonical point on the orbit; the polynomials 
which separate the generic orbits are the analogues of those used in the case n = 3. 

In the case of n = 2 this method can be used but there is some redundancy in the 
description we have given: the matrices Xn' = RijRi>j and Yjj> = RijRiji (using the notation 
of (0)) have the same eigenvalues and the matrix Rij is diagonal at the canonical point. 
In this case there are nine functionally independent invariants which specify the squares of 
the non-zero components of a , (5 and R at the canonical point on a generic orbit: trX n , 
ot T X m ~ x a and f3 T Y p ~ l (3 1 where n,m,p take the values 1,2,3. Additional invariants are 
needed to specify the signs of the non-zero components. The five invariants a.(Xa) A (X 2 a), 
{3.(Yj3)A(Y 2 j3) and aX r ~ 1 Rf3 ) r = 1, 2, 3, are sufficient to determine these signs for generic 
orbits and hence separate these orbits. In fact, using slightly different arguments, one can 
show that, in this case, one can reduce the number of polynomial invariants to ten, namely 
trX, trX 2 , detR, a T X T ~ X a, a T X r X R(3, r = 1,2,3 and Ag, which are subject to a single 
relation expressing A\ as a function of the other invariants. 

The general idea of investigating canonical points on orbits in the way we have described 
is also appropriate for higher spins, but the situation is somewhat more complicated. Con- 
sider the example of two particles of spin one in which case the unitary group under which 
p transforms is SU(3). p may be written as 

p = — 1 8 <g> 1 8 + ai T % ® 1 8 + fiiU ® Ti + R^ ® Tj (36) 

where Tj, i — 1...8 are representatives of a basis for the Lie algebra of SU(3) in the adjoint 
representation and Is is the 8x8 identity matrix. However, the adjoint representation of 
SU(3) is equivalent not to 5*0(8) but to an eight dimensional subgroup of it; this means 
that we cannot transform p so that RR T and R T R are diagonal so the canonical form is 
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rather more complicated than in the case of spin- 1/2 particles. 

In summary we have shown how to calculate the number of functionally independent 
parameters needed to determine whether or not two density matrices are locally equivalent. 
We have also shown how to characterise the generic classes of locally equivalent density 
matrices of n spin- 1/2 particles by two methods: (a) by finding an explicit set of non- 
polynomial invariants (the components of the density matrices at the canonical points on 
the orbits) and (b) by finding an explicit finite set of polynomial invariants. These methods 
work for generic density matrices; in a future publication we intend to give a systematic 
method for characterising classes of locally equivalent non-generic density matrices. In 
particular this will give a basis for the ring of invariants. We note that the canonical point 
on certain types of non-generic orbit has non-trivial stability group; this is a signature that 
density matrices on this orbit have special types of entanglement [0. 
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